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ABSTRACT 

Context. Binaries with hot massive components are strong X-ray sources. Besides the intrinsic X-ray emission of individual binary 
members originating in their winds, X-ray emission stems from the accretion on the compact companion or from wind collision. Since 
hot star winds are driven by the light absorption in the lines of heavier elements, wind acceleration is sensitive to the ionization state. 
Therefore, the over-ionization induced by external X-ray source strongly influences the winds of individual components. 

Aims. We studied the effect of external X-ray irradiation on hot star winds. 

Methods. We used our kinetic equilibrium (NLTE) wind models to estimate the influence of external X-ray ionization for different 
X-ray luminosities and source distances. The models are calculated for parameters typical of O stars. 

Results. The influence of X-rays is given by the X-ray luminosity, by the optical depth between a given point and the X-ray source, 
and by a distance to the X-ray source. Therefore, the results can be interpreted in the diagrams of X-ray luminosity vs. the optical 
depth parameter. X-rays are negligible in binaries with low X-ray luminosities or at large distances from the X-ray source. The 
influence of X-rays is stronger for higher X-ray luminosities and in closer proximity of the X-ray source. There is a forbidden area 
with high X-ray luminosities and low optical depth parameters, where the X-ray ionization leads to wind inhibition. There is excellent 
agreement between the positions of observed stars in these diagrams and our predictions. All wind-powered high-mass X-ray binary 
(HMXB) primaries lie outside the forbidden area. Many of them lie close to the border of the forbidden area, indicating that their 
X-ray luminosities are self-regulated. We discuss the implications of our work for other binary types. 

Conclusions. X-rays have a strong effect on the winds in binaries with hot components. The magnitude of the influence of X-rays can 
be estimated from the position of a star in the diagram of X-ray luminosity vs. the optical depth parameter. 
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1. Introduction 

Binary stars with hot massive components are well-known X-ray 
sources. A high-mass X-ray binary (HMXB) is a system consist¬ 
ing of a massive luminous hot star and a compact object (either 
a neutron star or a black hole). In these systems a fraction of 
the primary star matter, coming from either a hot star wind or a 
Roche-lobe ove rflow, is trapped in the gravitational well of th e 
compact object (iDavidson & Ostrikeri 1973tlLamers et al.l 1976) . 
These objects belong to the most powerful stellar X-ray sources 
with highest X-ray luminosities Lx ~ 10 37 - 10 38 erg s' 1 . 

Be/X-ray binaries also constitute a hot star donor and a 
compact companion, however, the matter is accreted from the 
circumstellar disk of the Be-star primary in these objects (see 
lReidl20fTl for a review). The X-ray luminosities of Be/X-ray 
binaries may reach up to Lx ~ 10 !7 erg s“'. Finally, the X- 
ray emission may originate in the w ind collision in binaries 
with non-degenerate hot comp onents (iPrilutskii & Usovlll976l: 
ICooke et alJfl978fclPittardll2009l) . This kind of interaction results 
in rel atively weaker X-ray sources with Lx ~ 10 32 - 10 33 erg s -1 
(e.g. JSana et al.ll2006bl : lAntokhin et alil20C)8l) . 

The hydrodynamics of the circumstellar matter in binary sys¬ 
tems is very complex and, in fact, constitutes a ti me-dependent 
problem of radiativ e hydrodynamics (iFriend & Caston [l982; 
Blon din et alJ 1 199ft iFeldmeier et al.l 1 19961) . Therefore. most 
studies concentrate on the dynamics of the matter leading to the 
X-ray generation. However, there is also a feedback effect of the 
produced X-rays on the stellar wind of hot component(s). 


The feedback eff ect is connected with X-ray photoionization 
of the stellar wind dMac Gregor & Vitellolll982l : 1van Foon et al.l 
1200 lb IWatanabe et ah 120061) . Because the stellar wind of hot 
stars is driven by light absorption in lines of heavier ele¬ 
ments (e.g., carbon, nitrogen, and i ron), photoionization may 
affect wind driving (IWatanabe et al.12006 ). This may even re¬ 
sult in wind stagnation and possibly fall back on the source star 
dKrticka et al .112013) . 

The problem of the X-ray photoionization should be treated 
using 3D hydrodynamical simulations coupled with radiative 
transfer and statistical equilibrium equations. These kinds of 
calculations are likely beyond the possibilities of current com¬ 
puters. To make the problem more tractable, one either uses 
hydrodynamical si mulations with a simp l ified form of the ra- 
diative force (e.g. IFriend & Castor! Il982b iBlondin et alJ 1 19901 : 
Feld meier et alJfl996b lHadrava & Cechurall2012l: iParkin & Siml 
20131) or ID stationary wind models that account for the in¬ 
fluence of X-ray photoionization on the level populations (e.g., 
IKrticka et ali [2012 1. Time dependent hydrodynamical simula¬ 
tions can explain the wind structure in binaries in detail, but they 
cannot predict for which binary parameters the effect of X-ray 
irradiation becomes important. 

We concentrate on a detailed solution of statistical equilib¬ 
rium equations and neglect the time-dependent phenomena. We 
provide a grid of stationary ID wind models with external X-ray 
irradiation for a wide grid of O star parameters. The X-ray irra¬ 
diation is parameterized by the corresponding X-ray luminosity 
and the distance of the X-ray source from the star with wind. 
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The wind is only calculated in the direction towards the X-ray 
source. This approach is sufficient to infer the effect of X-ray ir¬ 
radiation on the radiative force and yields a reliable estimate of 
the wind velocity in the direction of the X-ray source. 


one can introduce the ionization parameter 

£0) = J L v e ~ Tv<r) dv, (5) 


2. Influence of X-rays on the wind ionization 


The influence of the X-rays on the ionization state of the wind 
is stronger the lower the X-ray optical depth between the X-ray 
source and a given point in the wind. The radial optical depth 
between the X-ray source at radius D and point at radius r is 
given by 


r y (r) = 



K v (r')p(r') dr' 


( 1 ) 


with distances measured with respect to the centre of wind¬ 
losing star. The absolute value in this equation ensures that the 
optical depth is positive even for r > D. The X-ray optical 
depth T v (r) depends on density, which follows from the continu¬ 
ity equation as p{r) = M/(Anr 2 v r ), where M and ty are the wind 
mass-loss rate and the radial velocity, respectively. Assuming 
that the X-ray mass-absorption coefficient k v is spatially con¬ 
stant and that the wind already reached its terminal velocity u m , 
Eq. £0 can be written as 


Ty(r) 


KyM 

4-nVoo 


1 

D 


( 2 ) 


The wind opacity in the X-ray domain mostly stems from the 
direct and Auger ionization. Eq. ((2} may motivate one to plot 
stars in the Lx - r v diagram to es tima te the influence of X-rays 
on the wind ionization state (IKrticka et al.ll2012L Fig. 4). 

Inspired by Eq. (0 we introduce a non-dimensional X-ray 
optical depth parameter 

‘M a, 

Voc\R* Dj\ 10 8 M 0 year _1 / 


to characterize the radial optical depth between the stellar sur¬ 
face and the X-ray source. This parameter estimates the influ¬ 
ence of X-rays on the ionization state in the region of the critical 
point where the wind mass-loss rate is determined. The criti¬ 
cal point is defined as a point where the speed of the Abbott 
waves, the fastest waves in the wind, is equal to the wind ve¬ 
locity (lAbbottll 19801: lOwocki & Rvbickill 19841) . If the matter is 
highly overionized close to the critical point, then the wind is not 
accelerated radiatively since higher ions are not efficient wind 
drivers (e.g. IKrticka et al.1120121) . This leads to wind inhibition. 
Therefore, we expect a strong influence of X-rays in stars with 
large X-ray luminosity Lx and low X-ray optical depth param¬ 
eter fx, whereas we assume weak X-ray influence in stars with 
low Lx and high tx • For medium values of these parameters, the 
X-rays may just affect the wind velocity and lead possibly to 
wind stagnation. 

A description of the influence of X-ray irradiation on the 
wind state via the optical depth parameter tx introduced in 
Eq. (0 can only be used for stars with a similar ionization and 
density structure. For a more complete description of this prob¬ 
lem, one also has to account for the rr dependence of the re¬ 
combination that follows the X-ray ionization (n is the electron 
number density). Therefore, with the mean intensity of the radi¬ 
ation ,/ , x at the distance d = \D — r\ from the X-ray source (see 
also lKrticka et alJl2012i Eq. 4), 


J x = 


\6n 2 d 2 


-Mr) 


(4) 


which estimates the influence of X-ray radiation on the ioniza¬ 
tion equilibrium of the wind. Here L x is the X-ray luminosity 
per unit of frequency and the integration goes over the X-ray 
domain. In the optically thin limit this gives the ionization pa¬ 
rameter £, ~ Lx/(nd 2 ) introduced bv lTarter et al.l (119691 see also 
[Hatchett & McCravI 19771) . where Lx — f L x dv. 


3. Description of the X-ray irradiated wind models 

For our calcula tions, we used spherica lly symmetric stationary 
wind models of lKrticka & Kubatl (120101) . The radiative force was 
calculated using the solution of the comoving frame (CMF) ra¬ 
diative transfer equation with occupation numbers derived from 
the kinetic equilibrium (NLTE) equations. The model enables us 
to predict the wind structure, including the wind mass-loss rate 
and terminal velocity, consistently from global stellar parameters 
only. 

The wind ionization and excitation state was calculated from 
the NLTE equations. Ionic models for the NLTE calculations 
are based on the Opacity an d Iron Project results dSeaton et al.l 


1992: Hummer et al. 19931) and on the data described by 


Pauldrach et al.l (1200 ll) . Part of the ionic mod els was adopted 


from TLUS TY model atmosphere input files (lLanz & Hubenvl 
12003L 2007 ). The line radiative force was derived using the ra¬ 
diative flux calculated from the CMF radiative transfer equation 
dMihalas et al.lll97~5h with the actual occupation numbers calcu¬ 
lated from the NLTE equations. The line opacity data used in the 
line force calculation were extracted from the VALD database 
(Piskunov et al. l 19951 Kupka et al.1 1999 ). The occupation num¬ 
bers derived from the NLTE equations were also used to calcu¬ 
late the radiative cooling and heating terms (iKubat et~akli 19991) . 
The emergent surface flux (inner boundary condition) was taken 
from the H-He spherically symmetric NLTE model stellar at¬ 
mospheres of lKubat) ( 20031 and references therein). The result¬ 
ing hydrodynamical equations, that is, the continuity equation, 
equation of motion (with CMF radiative force), and the energy 
equation (with radiative cooling and heating), were solved iter¬ 
atively together with NLTE and radiative transfer equations to 
obtain the wind density, velocity, and temperature structure. 

Because the velocity field may become non-monotonic in 
the presence of the external irradiation, we do not use the CMF 
procedure of the radiative force calculation directly. Instead, as 
in IKrticka et al.l d2012h . we calculate the ratio of the CMF and 
Sobolev line force (see lKrticka & Kubad2010l ) for a model with¬ 
out external X-ray irradiation and use this ratio to correct the 
Sobolev line force in the models with external X-ray irradiation. 
The numerical test showed that this approach introduces very 
minor difference (below 1%) between models with CMF line 
force and models with scaled Sobolev line force. Moreover, we 
use the model atmosphere flux for the calculation of the Sobolev 
line force. Therefore, we neglect nonlocal radiative coupling 
betwe en absorption zones, which occurs in the non-monoto nic 
flows dRvbicki & Hummenll978tlFeldmeier & Nikuttall2006l) . 

The influence of the secondary component is only taken into 
account by inclusion of external X-ray irradiation. This is sup¬ 
posed to originate in the wind accretion on the compact compan¬ 
ion or in the wind-wind collision. The X-ray irradiation is mod¬ 
elled by an additional term in the mean intensity J v in the form of 
Eq. 0. The optical depth r y (r) along a radial ray in the direction 
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Table 1. Stellar parameters of the model grid. 


Table 2. Parameters of the velocity law fits. 



Model 

T c ff 

R, 

M 



[K] 

1^0 1 

[M e ] 

supergiants (I) 

300-1 

30000 

22.4 

28.8 


375-1 

37 500 

19.8 

48.3 


425-1 

42 500 

18.5 

70.3 

main 

300-5 

30000 

6.6 

12.9 

sequence(V) 

375-5 

37 500 

9.4 

26.8 


425-5 

42 500 

12.2 

45.0 


Model 

Vi [kms *] 

i>2 [kms *] 

i>3 [kms *] 

r 

300-1 

2280 

-680 

0 

-610 

375-1 

2700 

-980 

0 

-6300 

425-1 

2720 

0 

-380 

-18600 

300-5 

2680 

-130 

0 

-11700 

375-5 

3380 

-1320 

0 

-6900 

425-5 

2770 

0 

-510 

-13700 


of the X-ray source is given by Eq. ©. The frequency distribu¬ 
tion of emergent X-rays L x is for simplicity approximated by 
the power law L x ~ v _1 from 0.5 to 20 keV fcf. lWatanabe et al.l 
120061 ). The total X-ray luminosity L x = f L* dv and the distance 
of the X-ray source D are free parameters of the model grid. 

The absorption coefficient and the density in Eq. © can be 
derived directly from the actual model. However, to simplify the 
calculation of the external irradiation term Eq. © in the presence 
of the kink in the velocity law (see below), for the calculation of 
/ x we use density and absorption coefficient in the form of 

M 

P0) ~ 47rr 2 u(r)’ 
y(r) = min(D(r), uidnk), 

Kv(r) = 

where k x is the depth-independent approximation of X-ray opac¬ 
ity, fkink is the velocity of the kink (see below) if it is present 
in the models, and otherwise u k j nk = oo. The fits to the wind 
velocity v(r ) and absorption coefficient k x are derived from the 
model with no external irradiation. The wind velocity is fitted as 
(tKrticka & KubafeOl lb 


( 

r *\ i 

R *\ 2 ( 

R *\ 3 

H 

1 _ y) +l,2 ( 

+ "H 

) J 


x 


1 1 - exp 



(7) 


where Vi, t> 2 , '’ 3 , and y are free parameters of the fit. The X- 
ray opacity per unit of mass averaged for radii 1.5 R, - 5 R* is 
approximated as 


log 


1 cm 2 g _1 


min(fli log/f + b\,\ogao) A < A\, 

«2 log/l + bo, A > Ai, ''' 


where Ai = 20.18. The parameter A is non-dimensional, and has 
the same value as the wavelength in units of A. Here a lh a\, b \, 
« 2 , and bi are parameters of the fit. 

The models were calculated in the direction towards the X- 
ray source. The influence of the X-ray irradiation is strongest in 
this direction and it can be expected to decrease for the rays with 
increasing angular distance from the source. 


4. Calculated wind models 

The adopted grid of stellar parameters corresponds to O stars in 
the effective temperature range 30 000 - 42 500 K. The stellar 
masses and radii in Table Q] were calculated from the effective 
temperature using empirical relations of iMartins et al.l ( 2005a ). 
These are, together with the X-ray luminosity and the distance 
from the X-ray source, the parameters of the model grid. The 


Table 3. Parameters of the opacity fits. 


Model 

a 0 

a\ 

b\ 

a 2 

k>2 

300-1 

215 

2.619 

-0.917 

2.655 

-1.478 

375-1 

206 

2.558 

-0.857 

2.665 

-1.516 

425-1 

191 

2.482 

-0.794 

2.620 

-1.530 

300-5 

210 

2.618 

-0.917 

2.649 

-1.466 

375-5 

200 

2.516 

-0.819 

2.670 

-1.534 

425-5 

197 

2.497 

-0.804 

2.647 

-1.537 


X-ray luminosities and the X-ray source distances were selected 
to correspond to typical hot star X- ray sources. Wind param eters 
for the stars in Table©are given in lKrticka & Kubal (1201 2ll . The 
parameters of the velocity law fit Eq. © and opacity fits Eq. © 
are given in Tables [2] and [3] respectively. 

The resulting wind models are given in Figs. ©-© Here we 
plot the dependence of the radial wind velocity on radius in the 
wind models for individual studied model stars from TableQ]for 
different X-ray luminosities and X-ray source distances. 

For a low X-ray luminosity L x , or a large X-ray source 
distance D, the influence of the external X-ray source on the 
wind model is only marginal. In this case the X-ray ioniza¬ 
tion influe nces only trace ionization states (via direct and Auger 
ionization, MacFar lane et alJfl994t IPauldrach et al.ll 1994h . This 
modifies the emergent spectra because other trace ions with 
a high degree of ion ization appear (ISnow & Mortonl 119761 : 
ILamers & Snowl 1 1978 1, but the wind structure remains unaf¬ 
fected. Wind models for supergiant stars 375-1 and 425-1 with 
L x = 10 31 ergs -1 (see Figs. 0 and© are examples of models 
with negligible influence of X-ray irradiation. Because new ion¬ 
ization states appear that may drive the wind, the wind terminal 
velocity may be slightly higher dKrticka & Kubati2009l) . For ex¬ 
ample, this explains why 300-1 model with L x = 1 0 37 erg s _1 
and D = 5000 R e has higher terminal velocity than the same 
model, but with D = 7000 R a (see Fig.©. 

For higher X-ray luminosities L x (or a lower X-ray source 
distance D ), the effect of X-ray ionization becomes stronger. The 
ionization states that drive a wind in a case without X-ray irradi¬ 
ation become less populated close to the X-ray source, which 
results in a decrease of the radiative force and the wind ve¬ 
locity in the outer parts of the wind (e.g. 300-1 wind models 
with L x = 1 o 38 erg s ' 1 , Fig.©. The decrease of the radiative 
force may be so strong that the accelerating wind solution is 
not possible any more, and the wind velocity switches to de¬ 
celerating overloaded solutions (IFeldmeie r & Shlosmanl 12000: 
iFeldmeier et al.l 12008 ) with a typical kink in the velocity pro¬ 
file. This leads to a non-monotonic velocity law and decrease 
in the terminal velocity, which were deduced from observations 
of H MXBs bv iKaner et al.l (fl993h and calculated for Vela X- 
1 by iKrticka et al.l (12012lk The decelerating regions of the kink 
may be followed by accelerating regions if the kink is weak (e.g. 
300-1 wind models with L x = 10 32 ergs~' in Fig. ©or 425-1 
models with L x = 10 34 ergs~' and D > 25 R 0 in Fig.©. For 
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Fig. 1. Regions with different effect of the X-ray irradiation in the diagrams of X-ray luminosity (Lx) vs. the optical depth parameter 
(fx)- Graphs are plotted for individual model stars in Table [I] Individual symbols denote positions of: models with negligible 
influence of X-ray irradiation (black plus, +), models where X-ray irradiation leads to the decrease of the wind terminal velocity 
(red cross, x), non-degenerate components of HMXBs from Table[4](filled circles, #), and individual components of X-ray binaries 
from Table0(empty circles, O). The regions of the Lx - fx parameters that lead to the wind inhibition are denoted using the shaded 
area (H). 


stronger kinks, the switch back to the accelerating solution is 
missing (e.g. 375-5 models with Lx = 10 35 ergs~' in Fig. [6]). 
However, our calculations simplify the treatment of the radiative 
force in the kink because we do not account for the complex line 
resonances that occur in this case (e.g. lFeldmeier et aPl2008h . 


With decreasing source distance D (or increasing X-ray lu¬ 
minosity Lx) the kink position moves towards the stellar surface 
(see model 300-1 with Lx = 10 34 erg s~ 1 in Fig. [T|. If the kink 
occurs for velocities lower than the escape speed, the wind never 
leaves the star and the wind material falls back to the star (for ex¬ 
ample, this happens for the 300-1 model with Lx = 10 37 erg s~ 3 
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Fig. 2. Parameters of all models from Fig. Q] plotted in the diagram of X-ray luminosity vs. the ionization parameter. The ionization 
parameter Eq. (0) was evaluated at the wind critical point with radius r cr ; t . The influence of X-rays is denoted using different symbols. 
Blue plus signs + denote the weak influence of X-rays (possibly just on the ionization equilibrium), grey crosses X denote the 
strong influence of X-rays leading to the decrease of the terminal velocity, and empty red circles O denote models with wind 
inhibition. The individual symbols are slightly vertically shifted for a better readability. The extension of individual regions is also 
denoted using coloured areas. 


and D = 200 R e ). Our models only describe an accelerating part 
of the solution and the fall back has_to be studied using time- 
dependent models (see lPorter & Skouzal 19991 for a similar situ¬ 
ation). 

In the case when the kink reaches the critical point, the X- 
rays start to significantly influence the wind mass-loss rate. The 
wind may be completely inhibited by the X-rays. In some cases, 
it is still possible to find a stationary solution (with the kink now 
being a critical point, not shown here), however, with a mass-loss 
rate typically lower by one or two orders of magnitude. Time- 
dependent models are needed to study this problem. 

Derived results are summarized in the diagrams of Lx vs. tx 
in Fig. U] Each diagram indicates results for one of the model 
stars listed in Table 0 For a given star (with fixed M and (;«,), the 
optical depth parameter tx in Eq. (0) only depends on the X-ray 
source distance D, and tx is an increasing function of I). Wind 
models with a negligible influence of X-rays appear in the right 
bottom part of diagrams of Lx vs. tx (high tx and low Lx). Wind 
models in which the X-ray source causes a decrease in terminal 
velocity (denoted by the x symbol) appear in a diagonal strip of 
the diagram of Lx vs. tx- A forbidden area, where the X-rays 
lead to wind inhibition, is located in the upper left corner of the 
diagrams of Lx vs. tx diagrams. 

The influence of the X-rays may be simply described by 
the ionization parameter i) introduced in Eq. (01. This is shown 
in Fig. 0 where we plot the positions of all studied models 
in the diagrams of Lx vs. £. It follows that for low values 
of the ionization parameter, £ < 0.01 ergs -1 cm, the X-rays 
may influence the wind ionization state, but their influence on 


wind structure is negligible. The X-rays significantly influence 
the winds for £ > 0.01 ergs -1 cm. Large ionization parameters 
£ > 1 - 10ergs -1 cm lead to the wind inhibition. These param¬ 
eter values correspond to the forbidden area in the diagrams of 
Lx vs. tx- 

The borderline between the ionization parameters that lead 
to wind inhibition and those that do not is not the same for all 
model stars in Fig. 0 There are different reasons for this, the 
most important is likely the influence of the energy distribution 
of the irradiating X-rays. The X-rays become harder with in¬ 
creasing optical depth between the source and a given point (be¬ 
cause of frequency-dependent opacity, see Eq. 0. leading to a 
smaller influence of the X-ray irradiation. Therefore, a distant 
obscured source has a smaller effect on wind ionization than a 
closer and weaker unobscured source with the same £. The wind 
optical depth in model 300-5 is so low that there is no signifi¬ 
cant obscuration even for the distant sources. Consequently, the 
X-rays already influence the wind for £ ~ I erg s -1 cm. 

5. Comparison with observations 

5. 1. High-mass X-ray binaries 

The existence of the forbidden area in the diagram of Lx vs. tx 
can be tested against the observed parameters of HMXBs. For 
this purpose we searched the literature for available parameters 
of wind-powered Galactic HMXBs (these are listed in Table 0. 
In this table we restricted ourselves to HMXBs with primary ef¬ 
fective temperatures 7 d r > 25 000 K to avoid the bistability jump 
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Table 4. Parameters of HMXBs with neutron star or black hole companion 


Binary 

Sp. Type 

log (L/L G ) 

Teff [K] 

R[R a ] 

M[M e ] 

a [R 0 ] 

L x [erg s~‘] 

M [M 0 year 1 ] 

fx 

Reference 

2S 0114+650* 

Bllae 

5.61 

24000 

37 

16 

54 

1.1 x 10 36 

6.8 x 10“ 7 

0.9 

1, 2, 37 

Vela X-l ' 

BOIa 

5.63 

27000 

30 

23.5 

53.4 

3.5 x 10 36 

7.5 x 10- 7 

1.0 

3, 4,5 

4U 1700-377^ 

06.5Iaf 

6.02 

40200 

21.2 

58 

37 

2.2 x 10 36 

4.6 x 10~ 6 

4.6 

6, 7,8 

IGR J16418-4532 

08.51 

5.69 

32800 

21.7 

30 

31.6 

2 x 10 36 

9.8 x 10~ 7 

1.0 

9, 38 

IGR J18029-2016 

B1I 

5.59 

32500 

19.8 

20.2 

33.1 

2 x 10 36 

6.3 x 10- 7 

1.1 

10, 29 

IGR J16479-4514 

08.5Iab 

5.52 

31000° 

20 

35 

31.6 

1.09 x 10 34 

4.5 x 10- 7 

0.4 

11 

IGR J17252-3616 

BOIa 

5.79 

30000 

29 

15 

51 

1.6 x 10 37 

1.5 x 10~ 6 

11.3 

12, 13,40 

IGR J18483-0311 

B0.5Ia 

5.57 

24600 

33.8 

33 

95.7 

3.7 x 10 35 

5.7 x 10- 7 

0.8 

14, 15,42 

IGR J18450-0435 

091a 

5.58 

30000 

23 

30 

72 

7 x 10 35 

6.0 x 10- 7 

1.2 

16, 17, 43 

X Per' 

BOVe 

4.69 

32000° 

12° 

15.5 

420 

2.7 x 10 33 

2.2 x 10“ 8 

0.1 

18, 30 

IGR J11215-5952 

B0.5Ia 

5.73 

24700 

40 

29 

80* 

3 x 10 36 

1.2 x 10~ 6 

1.5 

19, 20 

PSRB1259-63 /c 

09.5Ve 

4.9 

34000 rf 

8.1 rf 

10 

1100 

3.5 x 10 34 

5.5 x 10- 8 

0.4 

21, 22 

IGR J19140+0951 

B0.5I 

5.47 

28000° 

23.2° 

25.4° 

62.5 

3 x 10 35 

3.6 x 10- 7 

0.7 

23 

4U 2206+54 

09.5Ve 

4.59 

30000 

7.3 

16 

53 

1.8 x 10 35 

1.4 x 10- 8 

0.1 

24 

4U 1907+09 

08.5Iab 

5.68 

29760 

26.2 

26.0 

54 

2 x 10 36 

9.6 x 10- 7 

1.6 

25 

LS 5039^ 

06.5V 

5.19 

37500 

9.3 

22.9 

34.5 

6 x 10 34 

2.0 x 10- 7 

0.7 

26, 39, 44, 45 

Cyg X-l/ 

09.7Iab 

5.57 

32000 

19.9° 

24.0 

42.4 

1.4 x 10 37 

5.7 x 10~ 7 

1.1 

27, 28, 41 

IGR J16465-4507 

09.51a 

5.69 

30000° 

26 

27.8 

124 

6.8 x 10 36 

9.9 x 10~ 7 

2.4 

31,32 

IGR J17544-2619 

091b 

5.53 

31000 

20.3 

26.5 

36.3 

1.7 x 10 35 

4.7 x 10~ 7 

0.7 

33,34 

IGR J16207-5129 

Blla 

5.4 

29000° 

20 

20 

40 

2 x 10 34 

2.6 x 10- 7 

0.5 

35,36 

XTE XI739-302 

08.5Iab 

5.71 

33000° 

21.9° 

30.8° 

173 

4.8 x 10 33 

1.1 x 10~ 6 

3.0 

46, 47 

IGR J17354-3255 

09.5Iab 

5.56 

30000° 

22.4° 

28.8° 

54 

2.6 x 10 36 

5.4 x 10- 7 

0.9 

48, 49 


Notes. Stellar parameters are taken from literature except for the mass-loss rate, for which we used fits of I Krticka & Kubiitj 120121) . ( °' Derived 
from the spectral type using the expressions of lMartins et alj (2005a). Periastron distance. (c) Possible disk accretion. <rf) Polar values.* 6 ' In the 
direction to the companion. Some alternative designations: 2S 0114+650 (V662 Cas), Vela X-l (GP Vel, HD 77581), 4U 1700-377 (V884 Sco, 
HD 153919), X Per (HR 1209), PSR B1259-63 (CPD-63°2495), LS 5039 (V479 Set), and Cyg X-l (V1357 Cyg. HD 226868). 


1 120081): (3) Ivan Kerkwiik et alJ 1 19951): (4) IStraizvs & Kurilienel 119811): (5)IWatanabe et all 
ikl 120021) : (8) Ivan der Meer et al.1 l2005h: (9) iDrave et alJ 120131) : ( TO) IMason et alJ 1201 lal) 


References. ( DlReie et al.l dl996l): (2) [Farrell et al 

ll2006[k (6) IColeiro & Chatvl l2013f): (7) Cla rk et al _ _ j _ _ _ _ 

(11) ISidoli et all 120131): 0 2) IMason et alj 1201 lbl); ( 13) iManousakis et alj 120121): (14) Romano et all 120101): (15) iRahoui & Chatvl120081) 
(16) Goossens et all 1 2013): (17) ISsu era et al. 1 20071): (18) Tomsick & Mutersnaugf l201 (]l): ( 19) Romano■ e t all (120 071): (20 ) iLorenzo et al 
1 2014 1) : (21) ISushch et al.l 20131): (22) ISierpowska-Bartosik & Bednarekl 120081): (23) IPrat et all 120081): (24)lRibo et all d2006l): (25)ICox et al 
1 2005|) : (26) ICasares et all 120051) : ( 27) iHerrero et al.l Il995h : (28 ) lHadrava & Cechural 120121): (29) iHill et al.l 120051): (30) iBerghbfer et a~ 

1 120141): (34)lPellizza et all 120061): (35) Tomsick et all 120091): 


(2008); (42) dearie et al.l <1200: 


( 19971): (31) IClark et alj 120101): (32) lLa Parola et alj l2010h: (33) IDrave et al ... 

(36)lBodaghee et al.ll201(T): (37)lHall et al.ll2000l): (38)lchatv et ahl(2008l) : (39)fHerrero et ahlll2002l) ; ( 40)lThomnson et alJl2007l) : (4~I) 


_ _ _ v .^-^jYar^emil 

(43) IZurita Herns & Walter! 120091): (44 ) iMoldon et alj 120121) : (45) iMcSwain et alj 120041) : (46) IDrave et al 


(2010); (47) iBodaehee et al J (hoi ll) : (48)Hg uera et al.l 1201 ll) : (49)ISg uer 1201 A . 


(iPauldrach & Pulslll990t IVink et al.lll999l) . which may compli¬ 
cate the estimate of wind parameters for cooler stars and de¬ 
serves special study. Our list of HMXBs with known parameters 
is likely incomplete, but it is a representative subset that charac¬ 
terizes the complete sample. 

Stellar and binary parameters (spectral type, 'L cl |, R, M, semi¬ 
major axis a, and Lx) in Table 0 were taken from the literature 
(however, see the notes below the table). We typically used the 
highest Lx in the case of variable X-ray sources. The mas s-loss 
rate was derived using Eq. (1) of [Krticka & Kubdtl (12012l) from 
the stellar luminosity L and the luminosity class. 

We included the parameters of HMXBs in the Lx - fx plots 
in Fig. [I] We assumed that the X-ray source distance D is equal 
to the semimajor axis a. The positions of all HMXBs in the di¬ 
agrams of Lx vs. fx lie outside the forbidden area. Moreover, 
many HMXBs appear just at the border between the forbidden 
area and the diagonal strip with a strong influence of X-rays on 
the wind. This indicates that the winds of these stars are in a 
self-regulated state. In the self-regulated state an increase in X- 
ray lumino sity cau ses wind inhibition and therefore a decrease 
in Lx (IKrticka et al.ll2012i) . whereas with a lower X-ray luminos¬ 
ity the influence of X-rays on the wind becomes weaker and the 
wind is faster, which leads to a stronger X-ray production. The 
combination of the two effects may keep the location of the star 


in the diagram of Lx vs. fx in the vicinity of the border line with 
forbidden solutions. 

Some objects (e.g. IGR J16479-4514, IGR J17252-3616, and 
CPD-63°2495) show strong X-ray variability, which would ap¬ 
pear as a vertical shift of their position in the Lx - fx plots. A 
horizontal shift in Lx-fx plots is introduced in eccentric systems 
as a result of the orbital motion (e.g. IGR J11215-5952, CPD- 
63°2495). This shows that ionization conditions in the wind may 
change with time. 

Parameters of most HMXBs lie in the area with a strong 
influence of X-rays on the wind structure (see Fig. |T]», leading 
to the decrease in wind terminal velocity with respect to non- 
irradiat ed winds. Within the classical Bondi-Hoyle-L yttleton 
picture dHovle & Lvttletonll 19411 iBondi & Hovlel[l944l) this af¬ 
fects the accretion rate and the X-ray luminosity. This effect can, 
for example, explain, why the wind velocity required to explain 
the X-ray luminosity of 4U 2206+54 is significantly lower than 
the terminal velocity expected for this type of star dRibo et al.l 
120061 ) . 

The existence of the forbidden area in the diagram of Lx 
vs. fx may be one of the reasons, why the X-ray luminosity of 
IGR J16479-4514 is by two orders of magnitude lower than that 
estimated from the Bondi-Hoyle accretion theory (iSidoli et al.l 
120131) . This star with L x = 1.09 x 10 34 ergs~' and fx = 0.4 
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lies very close to the forbidden area. Moreover, during flares it 
reaches peak luminosities of L x = 6.4 x 10 34 erg s 1 , which cor¬ 
responds to the boundary of the forbidden area. 

Stars exhibiting Roche-lobe overflow may appear in the for¬ 
bidden area, since forces other than the radiative force drives the 
flow. This is like l y the case of Cen X-3, which for parameters 
from lSuchv et a l. (2008) and lAtovan et alJ (12002) gives f x = 0.3; 
this is too low for its X-ray luminosity Lx = 5.4 x 10 37 erg s _1 
(ISuchv et al.i 12008 ). There are also observ ational indications 
of the Roche-lo be overflow in this system (ISuchv et al.l [2008: 
Nai k etalJ2011h . 


We also tested additional HMXBs (not listed in Table |4j 
with T e ff < 25 000 K. Their wind parameters may be influ¬ 
enced by the presence of the bistability jump, and consequently, 
we applied the mass-loss rate predictions of I Vink et al.l d2000t) 
for these stars. For IGR J00 370+6122, with parameters from 
iGonzalez-Galan et al.l d2014l) . we derive f x = 3.3, which is for 
Lx = 2.5 x 10 35 erg s 1 well below the forbidden region in Fig. |T] 
A similar res ult was derived for Swift J17 00.8-4139. With pa- 
rameters from lChakrabartv et al.l (i2002l) an d lMason et al.l (l2012h 
we derive the optical depth parameter f x = 1.5, which is for 
Lx - 3 x 10 36 erg s 1 again well below the threshold in Fig. [I] 
The object IGR J16318 -4848, with L x = 1.3 x 10 35 ergs 1 


( Filliatre & Chatvi 20041) and f x = 0.5 (with parameters from 
IChatv & Rahoui 2012), lies just at the border of the forbidden 


5.2. Binaries with non-degenerate components 

The comparison with observations in the case of massive X-ray 
binaries with non-degenerate components is complicated by the 
fact that the X-ray source distance can not be uniquely deter¬ 
mined from observations. To overcome this, we assume that the 
wind collisional front is in equilibrium, which means that the 
momenta deposited per unit of time and surface by the wind of 
both components are equal (e.g. !Antokhin et al.ll2004l) . 

p 1 v\=p 2 vl, (9) 

where the subscripts 1 and 2 denote the wind parameters of the 
primary and secondary, respectively. Here we study the wind 
at the intersection of both components. Denoting the radial dis¬ 
tances from the individual star centres Z>i and D 2 and using the 
continuity equation, we find 


by X-ray emission. Therefore, each binary appears in Fig.[T|and 
Table [5]twice, once for the primary and once for the secondary. 

Most of the binary member parameters studied lie in the al¬ 
lowed part of the L x - f x plots in Fig. |T| This shows that the 
winds of both components are not inhibited by X-ray emission 
and therefore X-ray emission in these objects may originate in 
wind-wind collisions. Many parameters correspond to the region 
with a strong influence of X-rays on the wind velocity struc¬ 
ture. These binaries may be in a self-regulated state, but of a 
different type than we identified in the HMXBs. In binaries with 
non-degenerate components, an increase in X-ray luminosity L x 
causes a decrease in the wind velocity and therefore a decrease in 
Lx dParkin & Simil20i 3i). Similarly, the decrease in L x leads to 
increase in the terminal velocity, which subsequently causes an 
increase in L x . This may be one of the effects that contribute to 
the observational conclusion that most binaries are not stronger 
X-ray sources than correspo nding single stars (e.g. ISana et al.l 
2006bl: lAntokhin et al]l2008l) . 

Some of the secondaries’ parameters lie in the forbidden re¬ 
gion in Fig. [I] (these are denoted by a superscript b in Table [5}. 
This indicates that their wind may be inhibited by the X-ray 
emission. The emission may originate in this case from the col¬ 
lision of the primary wind with the secondary star. It is also pos¬ 
sible that ID models are not applicable for modelling the wind 
in these systems, or that some parameters are not correct (e.g. 
mass-loss rate). Another possibility is that the basic stellar pa¬ 
rameters (luminosities, masses, and radii) need revision. 

For some secondaries, the X-ray source distance D 2 calcu¬ 
lated from Eq. GSl is lower than their radius R 2 (see a super¬ 
script a in Table [5]). For these stars, the primary wind may reach 
the secondary and the secondary wind may be completely in¬ 
hibited in the direction toward the primary. Parameters of these 
secondaries are not plotted in Fig. [7] Our models may also be 
oversimplified in this case. For example, we neglect the ra¬ 
diative force from the secondary acting on the primary wind 
dGavlev et alJI 19971) . which may avert the secondary wind inhi¬ 
bition. Moreover, a significant part of X-rays may also originate 


in shocks in supersonic winds similar to the case of single stars 


(Lucv & White 1980; Owocki et al. 19881;; 

Feldmeier et all19971; 

Babel & Montmerle 1997; ud-Doula et al. 

2014). 


6. Implications for other massive X-ray sources 

6.1. Be/X-ray binaries 


D\ _ Mm, oc 

L>\ M 2 v 2 'Oo 


GO) 


with a = D 1 + D 2 . Here we assumed that both winds reached 
their terminal velocities. 

Stellar and wind parameters of individual binaries are given 
in Table [5] The list is based on the known binaries included in 
the iBerghofer et al . ( 199 7;) sample supplemented by several bi¬ 
naries found in the literature. Stellar and binary parameters are 
taken from the literature. The mass-loss rates were derived us¬ 
ing Eq. (1) of iKrt icka (12014 corrected for actual radius) for 
main-sequence B stars and Eq. (1) of IKrticka & Kubiitl (120121) 
for other stars. For the terminal velocity we assumed the rela¬ 
tion between the escape velocity v esc and the terminal velocity 
Vo, - 2 .6 Uasc. This relation follows from both obs ervations and 
theory (lLamers et al.ll 19951 IKrticka & Kubatll2004l) . 

We included the parameters of X-ray binaries in the L x - f x 
plots in Fig. Q] If the X-ray emission originates from the wind- 
wind collision, then any of the winds should not be inhibited 


The X-ray emission in Be/X-ray binaries originates from the ac¬ 
cretion of the Be star disk material on the compact source (see 
lReigl20Til for a review). The origin of the disk is still a matter of 
debate. Some theories propose that the disk is fed by the wind. 
We used the list of Be/X-ray binaries given in IKrticka et al.l 
( 2015 ) to test whether the radiatively driven wind can survive 
the strong X-rays source in Be/X-ray binaries. In many Be/X-ray 
binaries the compact source is located at large distance from the 
Be star, and, consequently, it can not disrupt its wind. However, 
in some Be/X-ray binaries the compact object is located in such 
proximity to the Be star, that it disrupts its wind. This is another 
argument that disfavours the wind origin of Be star disks. 

6.2. X-ray binaries in the Magellanic Clouds 

X-ray binaries in the Magellanic Clouds typically show 
large X-ray luminosities of the order of 10 37 - 10 38 erg s _1 . 
Because of lower metallicity of LMC and SMC compared 
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Table 5. Parameters of binaries with non-degenerate components. 


Binary 

Sp. Type 

log (L/L 0 ) 

T e g 

[K] 

R 

[*o] 

M 

[M e ] 

M 

[Moyear^ 1 ] 

Voo 

[kmT 1 ] 

fx 

a 

[*o] 

[erg s~‘] 

Reference 

HR 8281 

06.5V 

5.7 

41200 

13.9 

40.2 

2.0 x lO -6 

2290 

3.6 

40 

2.3 x 10 82 

1, 2, 3 


09V" 

5.02 

34900 

8.9 

22.4 

9.7 x 10~ 8 

2400 





HD 152219 

09.5111 

5.07 

31900 

11.2 

26.2 

8.0 x 10^ 8 

2320 

0.2 

37.1 

8 x 10 3 ° 

4 


B1V" 

3.74 

21800 

5.2 

10.3 

3.4 x 10~ 10 

2250 





HD 47129 

081 

5.35 

33500 

14.1 

54. 

2.0 x 10~ 7 

2980 

0.4 

128.5 

8.34 x 10 32 

5, 30 


07.5111 

5.09 

33000 

10.8 

56. 

9.0 x 10- 8 

3560 

0.2 




HR 6187 

03V 

5.78 

44500 

13.14 

63 

2.9 x 10~ 6 

3090 

3.9 

38.0 

3.2 x 10 33 

1,6, 7,8 


05.5V" 

5.28 

39000 

9.54 

40. 

3.0 x 10- 7 

3100 





HR 6736 

03.5V 

5.75 

43850 

13.1" 

55 

2.6 x 10~ 6 

2860 

6.9 

4100 

1.2 x 10 33 

L9 


05V 

5.5 

40850 

11.2" 

36 

8.1 x 10~ 7 

2560 

2.8 




V712 Car 

031 

6.06 

43000 

19.3 

82.7 

5.4 x 10~ 6 

2720 

2.8 

53 

7.5 x 10 33 

10, 11 


031 

6.06 

43000 

19.3 

81.9 

5.4 x 10~ 6 

2700 

2.8 




HR 65 

09111 

5.2 

32000" 

13 

15 

1.5 x 10- 7 

1490 

0.3 

32 

7.2 x 10 32 

12, 13 


09111 

4.83 

32000" 

8.5 

21 

2.7 x 10“ 8 

2430 

0.0 




HR 5664 

07V 

5.27 

37500 

10.2 

21.6 

2.9 x 10- 7 

2090 

0.9 

33.8 

6.9 x 10 31 

1, 14 


09.5V 6 

4.64 

33000 

6.4 

12.4 

1.7 x 10- 8 

2140 

0.0 




HR 8406 

08.5111 

5.11 

32000 

11.7 

16.9 

9.7 x 10“ 8 

1750 

0.2 

25.5 

7.1 x 10 31 

1, 15 


09.5V" 

4.69 

28000 

9.4 

6.7 

2.2 x 10~ 8 

1230 





HD 93206A 

09.71 

5.7 

32000 

23 

40 

1.0 x 10~ 6 

1780 

2.0 

116 

3.4 x 10 32 

16, 29 


B2V" 

3.7 

20000 

5.9 

10 

2.3 x 10- 10 

2080 





HD 93206B 

08111 

5.3 

32600 

14 

14 

2.3 x 10~ 7 

1310 

0.6 

49 

1.2 x 10 32 

16, 29 


09V 

4.9 

32500 

8.9 

28 

5.6 x 10- 8 

2750 

0.1 




VI007 Sco 

07.5111 

5.4 

34000 

14.47 

27.24 

3.7 x 10~ 7 

1950 

0.5 

51.7 

4.4 x 10 32 

17, 18 


07111 

5.5 

34350 

15.99 

27.88 

5.9 x 10- 7 

1810 

0.9 




HR 1931 

09.5V 

4.64 

32600 

6.6" 

19 

1.8 x 10 -8 

2650 

0.1 

375 

2.1 x 10 32 

1, 19, 31 


BOV 6 

4.45 

30600" 

6.0" 

15.4 

2.6 x 10~ 9 

2520 

0.0 




HR 6535 

06V 

5.26 

39000" 

9.4 

32 

2.9 x 10~ 7 

2750 

0.8 

56.4 

5.2 x 10 32 

1, 20, 21, 22 


07V 

5.17 

37000" 

9.4 

32 

1.9 x 10- 7 

2790 

0.5 




HR 1899 

09111 

4.81 

32000" 

8.3 

23 

2.5 x 10~ 8 

2580 

0.1 

80 

7.4 x 10 32 

1, 21, 23, 24 


B0.5V 6 

4.19 

28000" 

5.3 

13 

1.8 x 10- 9 

2480 

0.0 




HR 1852 

09.511 

5.15 

31000" 

13 

11.2 

7.8 x 10~ 8 

1250 

0.3 

44 

6.5 x 10 32 

1, 21, 25, 26 


B0.5III 6 

4.2 

29000" 

5 

5.6 

1.5 x 10~ 9 

1640 

0.0 




V918 Sco 

07.51 

5.9 

35100 

24.3 

57 

2.7 x 10~ 6 

2000 

3.1 

91 

5.2 x 10 32 

1, 27 


09.71 

5.79 

30500 

28.1 

37 

1.5 x 10~ 6 

1430 

0.7 




HR 7767 

08.5111 

5.25 

32500 

13.4" 

23 

1.9 x 10- 7 

1900 

0.7 

620 

8.9 x 10 31 

1, 28 


B2.5V 6 

3.29 

20000 

3.7" 

9 

3.5 x 10~ H 

2500 

0.0 




DH Cep 

05.5V 

5.37 

44000 

8.31 

29.4 

4.5 x 10- 7 

2720 

0.9 

26.4 

5.4 x 10 33 

12, 32 


06.5V 6 

5.27 

43000 

7.76 

25.0 

2.9 x 10- 7 

2620 

0.5 




V1034 Sco 

09V 

4.78 

33200 

7.45 

16.8 

3.3 x 10~ 8 

2310 

0.1 

22.65 

3.2 x 10 31 

33,34 


B1V" 

3.87 

26200 

4.18 

9.4 

5.8 x 10 10 

2390 





VI294 Sco 

09V 

5.1 

34000 

10.3 

24.5 

1.4 x 10- 7 

2320 

0.4 

42.8 

6.9 x 10 31 

35 


09.7 V 

4.71 

31200 

7.8 

18.2 

2.4 x 10~ 8 

2370 

0.0 




HD 93205 

03V 

6.18 

49000 

17.2" 

45 

1.8 x 10“ 5 

1130 

60.8 

56.3 

3.4 x 10 33 

36 


08V" 

5.11 

36500 

9.0" 

20 

1.4 x 10~ 7 

2200 






Notes. Stellar parameters are taken from literature except for the mass-loss rate, for which we used fits of iKrticka & Kubdi <20l2h an d fKrt ic kit 
(1201411 . Alternative designations: HR 8281 (HD 206267), HD 152219 (V1292 Sco), HD 47129 (V640 Mon, HR 2422, and Plaskett star), HR 6187 
(HD 150136), HR 6736 (9 Sgr, HD 164794), V712 Car (WR 20a), HR 65 (AO Cas, HD 1337), HR 5664 (6 Cir, HD 135240), HR 8406 (14 Cep, 
LZ Cep, HD 209481), HD 93206 (QZ Car), V1007 Sco (HD 152248), HR 1931 (cr Ori, HD 37468), HR 6535 (V1036 Sco, HD 159176), HR 1899 
(t Ori A, HD 37043), HR 1852 (6 Ori, HD 36486), V918 Sco (HR 6164, HD 149404), HR 7767 (HD 193322), and HD 93205 (V560 Car). ( °> The 
primary wind reaches the secondary surface (Do < Ri )- (6) The wind parameters lie in the forbidden area. ^Derived from the effective temperature 
or spectral type using the expressions of lHarmaned <19881 ) for main-sequence B stars and Martins et all <2005al) for others. 


Refe rences. (1) iBeruhofer et aP <1997l): (2) Burkh older etaP d 19971): f3) ISticklandl dl995h : ( 4) ISana et al.l <2006al) : (5) 

(6) iMahvetal] <20121); (7) ISana etalj <20131): (8) ISanchez-Bermudez et aid <20131); (9) iRauw et all <20121): (10) 

(Tl) iMontes et alJ <2013l) : 0 2) IChlebowski et al] dl989l): 03) ISahade & Brandil d 19911): 04) iPennv et all <200ll): 05) iMahv et aid <201 ll): 


(16)jParkin etjdj 200); (17) Sana ct al 

(2003); (22) Chlebowski & OarmanvH 19911 


W (1 

!8)te 


Linder et al J ( 2006 ) 
Rauw et aH (2005); 


[2004 rO8)lMaveret all <20081) ; 0 9) ISimon-Dfaz et alJ d2011l); (20) iPachoulakisI dl996l); (21) IgI 
); (23) Ma rchenko et all <200(1): f24)lBagnuolo et alJ<200lh: (25)lHarvin et all<20021): (26)lMaver et al 

~) Li 


: 08}JtenBrummelaaret alJ <20111): (29)lLeung et all <19791): (30)lLinder et alJ <20081): (31)lNaiarro et al 

(32) ISturm & Simort <1994l) ~(33) ISana et al.1 <2005h : (34) lBouzid etaP <2005h : (35) ISana et all <2008h : f36) lAntokhina et all <200(1) . 


yer et al. 

J <201 lb ; 










































































J. Krticka et al.: The effect of the X-ray irradiation on the wind in binaries with hot component 


to our Galaxy, the wind mass-loss rates are expected to be 
also correspondingly lower (e.g. iBouret et alJ 1 2003 ). We tested 
three LMC and SMC binaries, LMC X-l, LMC X-4, and SMC 


parameters taken from literature (|Orosz et all |2009 ; 
lJ l2010t iRawls et alJl20llt iLi & Lill2014l) . In all these 


X-l, with 
iHung et al.l 

cases, the corresponding stellar parameters lie in a forbidden 
zone, indicating that the wind of these stars is inhibited in the 
direction toward the companion and the X-rays most likely orig¬ 
inate in the Roche lobe overflow. There is an observational sup¬ 
port for this indicatin g that many of these binar i es indeed l ikely 
fill their Roche lobes (Negueruela & Coell2002t IVal Baker et al.l 


2005;: lOrosz et~akll2009). 


6.3. X-rays from wind shocks 

The X-ray emission of single non-magnetic hot stars presum¬ 
ably originates in the sh ocks caused by the instability con¬ 
necte d with line driving dOwocki et allIl988k iFeldmeier et al.l 
1 19971) . There may be a feedback effect of the shock generated 
X-rays on the ionization structure of the wind. The typical X- 
ray luminosity scales with the stellar lu minosity. Lx ~ 10 -7 L 
dSana et alJl2006bt lAntokhin et al.ll2008h . however, simulations 
predict stochastic X-ray variability on some level. Taking as an 
example the predicted Lx variations of the 09.71b star C Ori A 
(in the order of JO 33 - 10 34 ergs -1 ) from the numerical simu¬ 
lations of IFeldmeier et al.l d 19971) and assuming the distance of 
X-ray source from the star 1.5 R, - 5 /?*, the derived value of 
the tx parameter is tx ~ 0.5 - 1.5. This lies in the region of 
a strong influence of X-rays on the wind ionization state (see 
Fig.QJfor a model 300-1). Consequently, shock generated X-rays 
may cause kinks in the velocity profile, which in turn may lead 
to moving absorption feat ures in the P Cygni line profiles (e.g. 
ICranmer & Owockilll996i) or increasing wind inhomogeneities. 

For some low-luminosity O stars the predicted mass-loss 
rates are si gnificantly higher than those inferred fro m observa¬ 
tions (e.g. feouret et al.l 120031 : iMartins et al.l l2005bl) . This dis¬ 
agreement is termed as a ’’weak wind problem”. This problem 
could possibly be explained as a cons equen ce of X-ray influ¬ 
ence on the mass-loss rate (Drew et aki119941) . For a typical lu¬ 
minosity of stars showing the weak wind problem, which have 
log(L/L 0 ) - 5, the X- ray luminosity can be inferred using the 
relation Lx ~ 10 7 L dSana et al.ll2006bt lAntokhin et al.ll2.0Q8h 
as log (Lx /1 ergs -1 ) ~ 31.6. If an effect of X-rays on the wind 
is to appear for these low X-ray luminosities, the X-ray source 
has to be in a very close proximity to the star (D < 1.5 R*, see 
Fig. |U). However, X-rays typically originate at larger distances 
IFeldmeier et al.1119971) . Moreover, in our calculations we as¬ 
sumed one spatially localized X-ray source, which is unlikely if 
the X-rays have to influence the entire wind. For a spatially dis¬ 
tributed X-ray source the effect of X-ray ionization is expected to 
be lower. Consequently, our results disfavour this explanation of 
’’weak wind problem”, and the weak wind problem seems to be 
more likely caused by an excessively large cooling length in the 


post-shock region (Lucv & White 1980: 

Martins et al.1 20051i 

Cohen et al. 2008;; Krticka & Kuball 2009; 

Lucy 2012j). 


6.4. Winds in active galactic nuclei 

The X-ray overionization may also be problematic for line- 
driven winds in the active galactic n uclei. Using typical parame- 
ters of these winds in simulations of iHigginbottom et al.l ( 2014 ) 
(Lx = 10 44 erg s and D = 10 16 cm), we derive the ionization 
parameter of the order of unity (neglecting absorption) even for 
the densest regions with p = 10~ 12 gcm -3 . Since typical wind 
densities are by two to four orders of magnitude lower, this sup¬ 
ports the conclusion that the winds of active galactic nuclei are 
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strongly inhibited, unless a strong absorption source is present 
jHigginbottom et ~aP l20i4i). 

7. Discussion and conclusions 

We study the effect of external X-ray irradiation on hot star 
winds in massive binaries. Hot star winds are driven by light 
absorption in the lines of heavier elements, therefore the wind 
acceleration is sensitive to the ionization stage of the wind. We 
used our own computer code to calculate NLTE wind models 
for typical O star and X-ray source parameters found in massive 
binaries to estimate the influence of external X-ray irradiation. 

The influence of X-rays is determined by the X-ray lumi¬ 
nosity, and the X-ray optical depth between a given point in the 
wind and the X-ray source, which follows from the distance to 
the X-ray source and X-ray opacity along the path. Therefore, 
our results can be interpreted in the diagrams of X-ray luminos¬ 
ity vs. the optical depth parameter (in this paper referred to as 
the diagrams Lx vs. tx)- The effects of X-rays are negligible in 
binaries with low X-ray luminosities or large distances between 
the binary component with wind and the X-ray source. In this 
kind of a case the external irradiation leads just to a modification 
of the wind ionization state, which has a relatively low influence 
on the radiative driving. 

With increasing X-ray luminosity, the influence of X-rays 
on the wind structure becomes stronger. There appears to be a 
typical kink in the velocity profile at the point where the X-ray 
overionization is so strong that it lowers the radiative force so 
much that is not able to accelerate the wind any more. With de¬ 
creasing distance between the star and the X-ray source the kink 
moves towards the star. If the kink approaches the critical point, 
where the mass-loss rates is determined, the wind becomes in¬ 
hibited by X-rays. Consequently, there is a forbidden area in the 
diagram of Lx vs. tx with high X-ray luminosities and low opti¬ 
cal depth parameters, where the X-ray ionization leads to wind 
inhibition. 

We compared the positions of known wind-accreting, high- 
mass X-ray binaries in the diagrams of X-ray luminosity vs. the 
optical depth parameter with respect to the location of the for¬ 
bidden area. All primaries of binaries with compact components 
lie outside the forbidden area, which is consistent with our wind 
models. Many of primaries lie close to the border of the for¬ 
bidden area, which means that the increase in X-ray luminosity 
would lead to wind inhibition. This indicates that the X-ray lu¬ 
minosities of HMXBs are self-regulated. 

The X-rays may also strongly affect the ionization state of 
circumstellar matter in colliding wind binaries. For many of 
them we predict a strong influence of external irradiation on the 
wind ionization state, pointing to wind self-regulation of another 
kind. In some cases the wind of secondary may be inhibited by 
the X-ray irradiation. We also discussed the application of our 
results for the X-rays generated by wind instabilities and for disk 
winds in active galactic nuclei. 

There are further observational tests of our models. The de¬ 
crease in wind terminal velocity may be detectable for individual 
stars as the shift of the blue edge of the UV P Cygni line profiles 
* Kjiper et~af) |2993f) or as the shift of the X-ray emission lines 
( Watanabe et alLoOhl) . Moreover, the kink in the velocity profile 
found in the models with strong X-ray ionization may lead to a 
localized stronger absorption in the blue part of the P Cygni line 
profiles resembling the discrete absorption c omponents (DACs, 
iKaner et afll 19941: Cranmer & O woe kill 19961) . 

Despite a good agreement between observations and theory, 
our models have limitations. First, we assume a spherical sym¬ 


metry and solve the wind equations only in the direction to¬ 
wards the X-ray source. While this may be applicable to esti¬ 
mate the overall influence of X-ray irradiation, we are missing 
many important observational effects that are connected with de¬ 
partures of the wind from s pherical symmetry (tFriend & Castor! 
119821: iBlondin et al.lll990l iFeldmeier et al.l 1 19961) . We also ne- 
glected the time-dependent phenomena, which may be important 
especially in binaries with eccentric orbits. Moreover, there is a 
mounting evidence that hot star winds show small-scale struc¬ 
tures (clumping). This may weaken the effect of X-ray ionization 
because of increased recombination (lOskinova et al.ll20T2h . 

Regardless of these limiting assumptions, our models clearly 
demonstrate the effect of X-rays on the ionization balance of 
stellar winds in massive binaries. Depending on the strength of 
the X-ray source and its distance from the wind losing-star, the 
X-rays are able, via their influence on the ionization structure, to 
change the wind driving force and, consequently, the properties 
of the wind. In a limiting case of a strong X-ray source, this 
leads to wind inhibition. This inhibition is supported by observed 
properties of binaries. 
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Fig. 3. Dependence of the radial velocity on the radius in the wind of 300-1 star model. Each graph corresponds to different lumi¬ 
nosity of additional X-ray source Lx labelled in the graph. Each radial velocity plot is labelled with a corresponding X-ray source 
distance D in R Q . 
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Fig. 4. Same as Fig.|3] except for a model star 300-5. 










































J. Krticka et al.: The effect of the X-ray irradiation on the wind in binaries with hot component, Online Material p 3 



r! ft * -1 



r/ ft* -1 



r/ ft* -1 



r/ ft* -1 



/-/ft* -1 



rl ft * -1 




Fig. 5. Same as Fig. [3] except for a model star 375-1. 
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Fig. 6. Same as Fig. [3] except for a model star 375-5. 
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Fig. 7. Same as Fig. [3] except for a model star 425-1. 
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Fig. 8. Same as Fig. [3] except for a model star 425-5. 
















































